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Abstract 

It is shown how the pre-exponential factor of the Feynman propagator for a large class of 
potentials can be computed using contour integrals. This is of direct relevance in the context 
of tunnelling processes in quantum theories. The prerequisites for this analysis are accessi- 
ble to advanced undergraduate students and involve only introductory courses in ordinary 
differential equations and complex variables. 
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I. INTRODUCTION 



The spectrum of certain differential operators encodes fundamental properties of 
different physical systems. Various functions of the spectrum, the so-called spectral 
functions, are needed to decode these properties. One of the most prominent spec- 
tral functions is the zeta function, which relates for example to partition sums, the 
heat-kernel and the functional determinant; see, e.g., Ref. 16. Zeta functions are often 
associated with suitable sequences of real numbers {^fc} feg ]N! which, for many applica- 
tions, are eigenvalues of Laplace-type operators. In generalization to the zeta function 
of Riemann, 



c*0O = E*~'> 

k=l 



one defines 



c(*) = EV, ( 2 ) 

fc=l 

where s is a complex parameter whose real part is assumed to be sufficiently large such 
as to make the series convergent. 

To indicate how zeta functions relate to other spectral functions, let us use the 
functional determinant as an example because it is going to be the focus of the article. 
For the purpose of relating zeta functions and determinants assume for the moment 
that we talk about a sequence of finitely many numbers {Xk}k=i- Considering them as 
eigenvalues of a matrix L, we have 

n 

det L = Y[ Afc, 
k=i 



which implies 

n d 
In det L = E In Afc = 



ds 



E K s - 

s=0 k=l 



In the notation of Eq. (T5]) this shows 

In det L = -C'(0) or det L = e~^°\ (3) 

When the finite dimensional matrix is replaced by a differential operator L having 
infinitely many eigenvalues, in general Yl^Li ^k will not be defined. However, as it 



turns out for many situations of relevance, definition ([3]) makes perfect sense and has 
found important applications in mathematics and physics; for the first appearance of 
definition see Refs. 10, 13, and 24. 

In recent years a contour integral method has been developed for the analysis of zeta 
functions,- 1 ^ which, although applicable in any dimension and to a variety of spectral 
functions, shows its full elegance and simplicity when applied in one dimension, and 
when applied to functional determinants. One of the main reasons for the relevance 
of determinants is the fact that the evaluation of the Feynman propagator involves 
this quantity. The probably most important field of application of functional deter- 
minants deals with tunnelling processes in quantum mechanics, quantum field theory 
and quantum statistics.— 1 ^ Because of its relevance a considerable number of articles 
in American Journal of Physics have been devoted to this topic; see, e.g., Refs. 1, 
5, 7, 14, and 15, and so we decided to also concentrate on this topic. Our aim is to 
show how and why a contour integral method is extremely well adapted for the eval- 
uation of in particular functional determinants. By demonstrating an additional way 
by which results may be obtained we enlarge the arsenal of techniques by a component 
already proven useful in recent research, see, e.g., Refs. 17 and 19. The probably 
most attractive feature of our approach is that all prerequisites are known to advanced 
undergraduate students of physics and mathematics. Namely, we only assume some 
working knowledge with Cauchy's residue theorem^ and some elementary facts about 
ordinary differential equations.— 

The outline of this article is as follows. We explain the basic ideas of our approach 
by looking at the zeta function of Riemann, and by evaluating Cj?(0)- This is identical 
to the evaluation of the functional determinant of a free particle in an interval with 
Dirichlet boundary conditions at the endpoints. We then will consider the case of 
particles in a harmonic oscillator potential previously considered in Refs. 1, 7, and 
15. Results will be trivially rederived. Finally, we show how particles moving in 
any potential (satisfying reasonable conditions) and obeying quite general boundary 
conditions can be analyzed. The Conclusions highlight the most important points of 
our contribution. 
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II. FUNCTIONAL DETERMINANT OF A FREE PARTICLE IN AN INTER- 
VAL 



A free particle in an interval is described by the operator j& together with some 
boundary condition. In the context of the Feynman propagator Dirichlet boundary 
conditions are quite common^ and this is what we first concentrate on. It will be 
convenient to make a rotation in the complex t-plane and to define t = —it. The 
resulting operator 

p. * 



dr 2 

in terms of r has positive eigenvalues. This is the relevant setting in the context of 
quantum tunnelling,-&22£L22i2£ cf. the last part of Section HVl 

So in order to evaluate the functional determinant associated with this situation we 
consider the eigenvalue problem 



4>ti{t) = A n n (r), 0„(O) = (j) n (L) = 0. 



dr 2 

The eigenf unctions have the form 

n (r) = asin(^Jx^r) + 6cos(y / A^r). 

The appearance of the cosine is excluded by the boundary value <f) n (0) = 0. The 
eigenvalues are then found from the equation 



sm(^X n L) = 0. (4) 
This condition is simple enough to be solved for analytically and one determines 

j / T17T \ ^ 

n (r) = asin( v /A n r), A n = ( — J , n G IN, 

with some normalization constant a. 

Although in this particular case it is of course convenient to have an explicit expres- 
sion for the eigenvalues, let us pretend the best we are able to obtain is an equation 
like (jl]), namely eigenvalues are determined as the zeroes of some function F(X). As 
we will see, actually this is as convenient as having explicit eigenvalues, but of much 
larger applicability. 



For the given setting the natural choice F(X) = sin(VAX) has to be modified as 
A = satisfies F(0) = 0. In order to avoid F(X) having more zeroes than there are 
actual eigenvalues we therefore define 

Note that F(X) is an entire function of A. The next step in the contour integral 
formalism is to rewrite the zeta function using Cauchy's integral formula. Given that 
F(X) = defines the eigenvalues A n , then 

has poles exactly at those eigenvalues. Furthermore, expanding about A = A n we see 
for F'(X n ) ^ that 

F'(X) F'(X - A n + A n ) F'(X n ) + (A - X n )F"(X n ) + ... J_ 

F(X) F(X - X n + A n ) (A - X n )F'{X) + (A - A n ) 2 F"(A n ) + ... A - A„ " 

and the residue at all eigenvalues is 1. (A variation of this argument shows that if m n 
is the multiplicity of X n , the residue of F > (X)/F(X) at X n is m n .) This shows, noticing 
the appropriate behavior of F(X) at infinity, that for 5Js > |, 

CpW = ^/«A-A™ (6) 

7 

where the contour 7 is shown in Figure 1. 

As is typical for complex analysis, the next step in the evaluation of a line integral is 
a suitable deformation of the contour. Roughly speaking, deformations are allowed as 
long as one does not cross over poles or branch cuts of the integrand. For the integrand 
in (151) , the poles are on the positive real axis, and there is a branch cut of A~ s which we 
define to be on the negative real axis, as is customary. So as long as the bahavior at 
infinity is appropriate, we are allowed to deform the contour to the one given in Figure 
2. 

In order to better see the | A | — > 00 behavior of F(X), let us rewrite the sine in terms 
of exponentials. We then have 

F(X) = — ^= (e l ^ L - e~ l ^ L 



A-plane 




FIG. 1: Contour 7 




A-plane 



FIG. 2: Contour 7 after deformation 



and for > | all deformations are indeed allowed. We next want to shrink the 
contour to the negative real axis as shown in Figure [3j 

As A approaches the negative real axis from above, A~ s picks up a phase (e t7T )~ s = 
e~ l7TS , whereas the limit from below produces (e~ t7T )~ s = e l7TS . Given the opposite 
direction of the contour above and below the negative real axis, contributions add up 
to produce a sin(7rs). Taking care of the same kind of argumentation in F(X) one 
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A-plane 



A = e^x 



A = e~ i7T x 



FIG. 3: Contour 7 

obtains 

Cp(s) = ^ Jdx x-±ln [1 - e-^ll • (7) 

Notice, that by shrinking the contour to the negative real axis a new condition for the 
integral to be well defined, namely 9Rs < 1, has become necessary due to the behavior 
about x = 0. 

Let us rest for a moment to stress the nice features of equation (jjj) for the evaluation 
of determinants. If the integral were finite at s = 0, an evaluation of the determinant 
would be trivial. From 

V 7T / «=0 l<£s«=°./ fife V 2 V^ L V/ 

it would amount to finding ln(...) at the limits of integration; even no integration needed 
to be done explicitly. Whereas this is exactly what occurs when considering ratios of 
determinants, see Section 4, for absolute determinants the situation is slightly more 
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complicated. The reason is that (|7j) is only well defined for | < 3?s < 1 and a little 
more effort is needed. Notice, the problem is caused by the x — > oo behavior which 
enforces the condition | < 3?s. To analyze this further we therefore split the integral 
according to Jq 1 dx + dx. Whereas from the above remarks it follows that Jq 1 dx can 
be considered to be in final form, the dx needs further manipulation. The pieces 
needing extra attention are 



I dx x 8 4-hxe^ L = — [ dx 
J dx 2 J 



x 



2s 



dx x s — In f — = 

dx \ 2Jx 



2 J ^ x x ~ 
i 



1 

'2s' 



This shows 

Cp(s) 



Lsimrs sin7rs sin7rs 



(2s - 1)tt 2stt 
i 



[ dx x s — In fl 

7T J dx V 



7T 



In 



\2y/x 



-2y/xL 



a form perfectly suited for the evaluation of Cp(0)- We find 



Cp(0) 



- In (l 



-2L 



+ In 



-2L 



-InL = -ln(2L). 



This, of course, agrees with the answer found from the well known values Cr(0) = ~\i 
&(0) = -iln(27r): 



CpW = E 



n=l 



L J 



Cp(2s) 



Cp(0) = 2 In ('£) Cfl(O) + 2&(0) = - In Q - 1h(2tt) = - ln(2L). (9) 



III. FUNCTIONAL DETERMINANT FOR PARTICLES IN A HARMONIC 
OSCILLATOR POTENTIAL 



Let u be the frequency of the harmonic oscillator, then the relevant operator to be 
considered is 



ho 



dr 2 



CO 



with Dirichlet boundary conditions imposed at the endpoints r = and r = L. Eigen- 
values are then determined by the implicit equation 



sm(sj\ n - uj 2 L) = 0. (10) 

Instead of looking at the determinant of Pho itself, let us now consider the ratio 
det (Pho) I det(P), where as before P = — J^, that is we consider the difference of 
the associated zeta function. Using the same strategy as before in Section 2, we have 

. , . ,. , . 1 f . d . / sin(VA - u 2 L) \/X \ 
0>,M - frW = 2tH I dX [ L (V /AL) \ zx-J>) ' 

the contour 7 still given by Figure 1. Deforming as before, we obtain 

Cft ( s ) - Cp(s) = / cb 2 — In , . 7 - v , (11) 

where sin(iy) = z sinhy has been used. The technically simplifying consequence of con- 
sidering ratios gets now apparent: as x tends to infinity, the behavior of the integrand 
has improved. In detail we have as x — > 00 

sinh(v^T^L) Vx LfVOT-^ 1 ~ e~ 2L ^ _ 1 u 2 L 

sinh( x /i^) V^T^ 2 v^T^ 2 l-e- 2L v^ 2^ 

and the integrand behaves like x~ s ~ 3 ^ 2 . Noting that the x — > behavior up to a 
proportionality constant is as before, we see that ( TTTT) is well defined for — | < < 1, 
in particular, it is well defined at s — 0. Thus, trivially, following along the lines leading 

to §]), 

CK o (0)-Cp(0) = -1u 



/ sinh ujL 



V J 

or, switching back to real time, replacing L = i(tf — ti), 

^ det P ho _ ^ sinh(iu;(t / - tj)) _ ^ sm(u(t f - tj)) 
detP iu(tf — ti) u(tf — ti) 

the well known answer; see, e.g., Refs. 11 and 15. 

Other boundary conditions can be dealt with basically with no extra effort. For 

example let us consider quasi-periodic boundary conditions as they have been analyzed 

for anyon-like oscillators.- 1 ^ In this case the boundary condition reads 

ML) = e te 0n(O), <f>' n (L) = e*V;(0), 



with 9 some real parameter; 6 = corresponds to periodic boundary conditions, 
whereas 9 = tt gives antiperiodic boundary conditions typical for fermions. The general 
form of eigenfunctions is 

(f) n (r) = asm (j X n - uj 2 r^j +bcos (j X n - iv 2 rj . 



The boundary condition produces the equations, use /i n = a/A„ — oo 2 , 

asin(/x n L) + 6cos(/i n L) = e^b, — fi n b sm(fi n L) + fi n a cos(fi n L) = e %e [i n a. 

Under the assumptions that /i n 7^ 0, which excludes periodic boundary conditions, this 
system represents the matrix equation 



0. 



sin(/i n L) cos(/i n L) — e %e \ 1 
cos(yU n L) - e ld - sm(jj, n L) J y 

This has a nontrivial solution if and only if the determinant of the matrix is zero, which 
after some simple manipulations gives the condition for eigenvalues as 

cos(/x n L) — cos# = 0. 

Following the steps of the previous calculation, denoting the operator with quasi- 
periodic boundary conditions as and P qp , the answer can essentially be simply 
i*6£icl off 

c^(o)-cr(o)--h( cosh 1 ( " L) - cose '), 

ho \ 1 — cos a j 

and agrees with Ref. 3. So in real time, 

detP^ ^ cos(k;(t/ — ti)) — cos 6 
det Pip~ 1 - cos 9 

For periodic boundary conditions an eigenfunction with zero eigenvalue occurs, namely 

the constant, and we comment on this situation in the conclusions. 



IV. FUNCTIONAL DETERMINANTS OF PARTICLES IN GENERAL PO- 
TENTIALS 

As the previous section made clear, the answer was obtained without ever worrying 
what the actual eigenvalues of the operator in question might be. The only information 
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that entered was the implicit eigenvalue equation ffTUj) . Is there any way an equation like 
([TO]) can be obtained for general potentials, such that the evaluation of determinants 
is similarly trivial as the previous one? The answer is yes and elementary knowledge 
of ordinary differential equations is all that is needed.— 

So let us say we were interested in the ratio of determinants of operators of the type 



2 



Pi = -^ + i2,(T), J = 1,2, 

where for convenience again Dirichlet conditions are considered. In the previous sec- 
tions Riij) = was chosen, but no additional complication arises for this more general 
case. Such ratios arise, for example, in the evaluation of decay probabilities in the 
theory of quantum tunnelling.-"^!^ Recall that if a quantum particle moves in a po- 
tential V(x) for which classically a particle is at rest at x — 0, and if x denotes the, 
say only, stationary point of the Euclidean action, then to leading order in h the decay 
probability per unit time of the unstable state is a multiple of, see Eq. (2.25) of Ref. 
9. 



det (~ + V"(m) 



-1/2 



Our contour integration method can easily handle such ratios. Indeed, as suggested 
by the previous examples, in order to evaluate det Pi/ det P2 in the general case, the 
contour integral to be written down should involve solutions to the equation 

p Aa(7) = ^,a(t), 

where A, for now, is an arbitrary complex parameter. As is well known, for continuous 
potentials Rj{r) there will be two linearly independent solutions and every initial value 
problem (j)j t \(0) = a, <//• A (0) = b, will have a unique solution. A contact with the original 
boundary value problem is established by imposing </>j,a(0) = 0; the condition on the 
derivative is merely a normalization and for convenience we choose (p'j A (0) = 1. The 
eigenvalues for the boundary value problem are then discovered by imposing 

<M£) = o, (13) 
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considered as a function of A. To see a little better how this works, consider the case 
i?2 = 0. The unique solution of the initial value problem described is 

Mt) = —7>r- 

Eigenvalues follow precisely from the condition 

ML) = 0. 

But having the implicit eigenvalue equation ( 1131) at our disposal, the calculation of the 
determinant is basically done! Arguing as below (jSJ) we write 

1 l„ >-. i 



C«W-Cp i M = ^-/^a--.„ 



7 



dX <p 2 ,\{L) 



7T J dX 02,-xW 

valid about s = because the leading behavior as x — > oo of <pj- x {L) does not depend 
on the potential Rj(r); as evidence see the analysis in Section 3. So as around (jSJ) 

and we obtain the Gel'fand-Yaglom formula^ 

detfi = 0!, O (L) 
detP 2 2i o (L)' 

The ratio of determinants is determined by the boundary value of the solutions to the 
homogeneous initial value problem 
/ d 2 



dr 



+ i? J (r)U,- (r) = 0, <M0)=0, 0;. O (O) 



Even if no analytical knowledge about the boundary value might be available, they can 
easily be determined numerically. 



V. CONCLUSIONS 



The main aim of this contribution was to show that the analysis of functional deter- 
minants for a large class of operators is accessible to advanced undergraduate students. 
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The only prerequisites are elementary ordinary differential equation theory and a basic 
course in complex variables. The beauty of the approach is that it is easily adapted to 
different situations. We have indicated how other boundary conditions than Dirichlet 
ones can be dealt with. Indeed, general boundary conditions can be considered along 
the same lines and generalizations of the Gel'fand-Yaglom formula can be obtained.— 

We have mentioned that the presence of zero eigenvalues adds some extra com- 
plication. The reason is that when deforming the contour to the negative real axis, 
some contribution from the origin may result. But again, a minor modification of the 
procedure allows for a complete analysis.— 

Even systems of differential equations can be considered with about the same 
effort.— 

An example where all of the above generalizations need to be considered is the study 
of transition rates between metastables states in superconducting rings. For this case, 
a 2 x 2-system with twisted boundary conditions needs to be analyzed; see, e.g., Refs. 
18 and 26. 

The many advantages of this approach described show that it is optimally adapted 
to the evaluation of determinants. Instead of struggling with the needed mathematical 
manipulations the students should be able to easily get a grasp of this technique and 
to concentrate on the underlying physics. 
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